Higher Order Collins Modulations in Transversely Polarized Quark Fragmentation 
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The Collins effect describes the modulation of the hadron production by a transversely polarized 
quark with the sine of the polar angle, between the produced hadron's transverse momentum 
and the quark spin. We employ a quark-jet model to describe multiple hadron emissions by such 
a quark, taking the Collins effect into account. The resulting hadron distributions exhibit modu- 
lation up to fourth order in sin((^) when only two hadron emissions are allowed, rising with any 
further increase in the number of emitted hadrons. These new effects are a direct consequence of 
the quark-jet mechanism for quark hadronization, which do not depend on the details of the model 
used for elementary hadron emission. The size and the sign of the higher order terms are directly 
connected with the probabilities of quark spin flip in the elementary emission process, with opposing 
sign favored and unfavored Collins functions only being generated if quark spin flip is preferential. 
Experimental studies of these effects should therefore provide a critical test of the quark hadroniza- 
tion mechanism, which in turn will lead to a deeper understanding of the transverse spin structure 
of hadrons. 
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Quark hadronization remains one of the most chal- 
lenging problems in hadronic physics, as it involves both 
long and short-range parton interactions. This process is 
quantitatively described by various fragmentation func- 
tions, which have a probabilistic interpretation. For ex- 
ample, the unpolarized fragmentation function, d\^^[z)^ 
can be interpreted as the probability density for an un- 
polarized quark q to emit a hadron h carrying light-cone 
momentum fraction z. These functions enter the expres- 
sions for the cross-sections of various hadronic processes 
in the factorization regime, such as semi-inclusive deep 
inelastic scattering (SIDIS) and e+e~ and hadron-hadron 
collisions, convoluted where necessary with the relevant 
parton distribution functions. Thus, a reliable knowledge 
of the fragmentation functions is crucial for disentangling 
the underlying parton distribution functions. 

In recent years there has been particular interest in 
studying scattering processes with transversely polarized 
probes and targets. Here the fragmentation of the trans- 
versely polarized quarks to spin zero particles, like pions 
and kaons, can be described by two functions: the famil- 
iar unpolarized fragmentation function and the naively 
time-reversal-odd Collins function [U |2] . 

The probability of the transversely polarized quark q 
to emit a spin-zero hadron h with transverse momentum 
Pj_, depicted schematically in Fig. [ij can be expressed 
as [3] 

D^,^,{z,Pl,^) = D\'\z,Pl) (1) 

ZTTlfi 

where rrih is the mass of the produced hadron. The trans- 



verse momentum dependent (TMD) unpolarized frag- 
mentation function is denoted by d(^^^(z,PJ), while 
H^^^^{z^P'j_) is the Collins function. Thus the Collins 
function modulates the distribution of produced hadrons 
with sin((/?), a clear signature that allows one to extract 
it from the experimental data. The experimentally mea- 
sured Collins functions from HERMES, COMPASS and 
JLab strongly suggest that the 1/2 moments 

of the unfavored Collins functions have a similar size and 
opposite sign to those for the favored ones for pions [4]- 
8 . We note that the fragmentation functions are called 
favored, where the produced hadron has a valence quark 
of the same flavor as the initial fragmenting quark, and 
unfavored otherwise. There has been extensive model- 
ing of the Collins function in the past [9-12 . Most no- 
tably, the spectator model calculations of Refs. [9HTT] 
provide a microscopic description of the fragmentation 
process. However, these calculations have been limited 
to the single hadron emission approximation, with the 
consequence that the unfavored fragmentation functions 
are zero in this approach. In addition, the description of 
the small z region for the favored fragmentation functions 
requires the introduction of a phenomenological quark- 
meson form factor. 

Recently, we employed the NJL-jet model of Refs. [13]- 
il6 to calculate the Collins fragmentation function within 
the quark-jet hadronization picture [17], using Monte 
Carlo (MC) simulations. In this picture, the initial quark 
undergoes a decay chain of hadron emissions, schemati- 
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cally depicted in Fig. |2j where a microscopic description 
of each elementary emission is provided by the effective 
quark model of Nambu and Jona-Lasinio (NJL) [18] . 

In this Letter we show that within this picture the 
fragmentation function of a transversely polarized quark 
acquires Collins modulations of higher order in sin((/p) as 
a direct consequence of the quark-jet hadronization pro- 
cess. These new terms are essential to describe the un- 
favored fragmentation functions. It appears that, even 
with just two hadron emissions in the decay chain, one 
requires a fourth order polynomial in sin ((/:?) to fully de- 
scribe the produced number densities. We stress that 
this effect is independent of the particular form of the 
model used to describe the Collins effect at the elemen- 
tary emission vertex but is a direct consequence of the 
multiple hadron emission mechanism. 

In the NJL-jet model of Matevosyan et al [T7| we used 
the tree level diagrams to calculate the unpolarized frag- 
mentation functions (the first term in Eq. ([T])), while for 
the Collins functions we used the spectator model of Bac- 
chetta et al. [10 . In that work, the interference term 
between a tree level amplitude and the amplitude with 
gauge link coupling was used to produce non-zero Collins 
function. 

A key element of this calculation was the evaluation of 
the quark spin-flip probability in each hadron emission 
step. We used light-cone spinors to describe the initial 
and final polarized quark wave- functions, and found that 
the quark spin non-flip and spin flip probabilities are pro- 
portional to 

|ai|2 ~ ll, |a2|2 r^ll + (M2 - (1 - z)M^f, (3) 

where Mi (2) is the constituent mass of the fragmenting 
(remnant) quark and is the transverse momentum of 
the remnant quark. (Note that this momentum cancels 




FIG. 1. Illustration of transversely polarized quark fragmen- 
tation. The quark momentum, k, defines the z-axis, with its 
transverse polarization spin vector, Sq, along the x axis. The 
emitted hadron has momentum P with component P± trans- 
verse with respect to the z-axis. The polar angle of P± with 
respect to Sq is denoted (p. 




FIG. 2. NJL-jet model including transverse momentum and 
quark polarization transfer. Here the (orange) double-lined 
arrows indicate the quark's spin. 

the transverse momentum of the emitted hadron with re- 
spect to the momentum of the fragmenting quark.) Here 
we clearly see that the quark spin-flip probability should 
be on average larger than the spin non-flip probability. 
This will be shown later to be crucial in describing the 
Collins function. 

We use the elementary fragmentation function for a 
transversely polarized quark calculated within the NJL 
model, along with the quark spin flip probabilities ex- 
plained above, as the input to a Monte Carlo simulation 
of the quark hadronization process within the quark-jet 
model - see Ref. [17] for a more detailed discussion. We 
used the and (p dependence of the resulting hadron 

number densities to extract the unpolarized and Collins 
functions for pions and kaons produced by light and 
strange quarks. The transverse momentum integrated 
number densities were perfectly described by a first or- 
der polynomial in sin((/p) and hence the constant term was 
interpreted as the integrated unpolarized fragmentation 
function, while the coefficient of the linear term as twice 
the 1/2 moment of the Collins function (c.f. Eq. ([2|). As 
perhaps the most remarkable result, we note that the un- 
favored Collins function was found to have the opposite 
sign and a magnitude comparable to the favored one for 
pions, as can be seen from the plot in Fig. |3] 




z 

FIG. 3. Fitted values for 2H^^^^^^ for pions produced by a 
quark from MC simulations with 6 emitted hadrons. 
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FIG. 4. Histogram of the values of Xdof for cill fragmentation 
functions of u quark fitted with polynomials of different order 
in sm((p) for MC simulations with 2 emitted hadrons. 



Even more interesting effects were observed when we 
considered the extraction of the TMD Collins function 
from the unintegrated number densities. We used a poly- 
nomial form in sin((/:?) to minimise ^ the (p 
dependence for fixed values of z and P^. Here, even 
for only two hadron emissions in the jet, the fits with 
just a linear form were insufficient. On the other hand, 
a fourth order polynomial in sin((/:?) provides an excel- 
lent fit everywhere. The plots in Fig. |4] illustrate the 
per degree of freedom (Xdof) ^^^^ polynomials of 
different orders for all the fragmentation functions of u 
quark (to seven pions and kaons in total) and for all the 
discrete values of z and in our simulations - a total 
of 7 X 100 X 100 = 7 • 10"^ data points. Clearly the fourth 
order polynomial yields Xdof vicinity of 1 for all 
the fits. A similar analysis for the integrated num- 
ber densities shows that a simple first order polynomial 
describes all the data perfectly well, meaning that this 
effect is only present in the TMD fragmentations. 

These higher order modulations can be attributed to 
two effects. First, in the process of multiple hadron emis- 
sions, the distribution of the n-th produced hadron is 
multiplied by an additional factor of (a + bsin{(p))^'^~^\ 
when averaged over all Monte Carlo simulations. This 
can be easily deduced from Eq. ([T]), as with each hadron 
emission the remnant quark number densities, relative to 
the fragmenting quark ((i^), are related to the number 
density of the emitted hadron (dq) by the momentum 
and flavor conservation law d^{l — z^ —p±) d^{z^p±_)^ 
where is the transverse momentum of the produced 
hadron with respect to the fragmenting quark. Thus 
these quark densities also get modulated with a first order 
polynomial form. The next emitted hadron's distribution 
is modulated not only by the elementary Collins function, 
but also by the number density of the fragmenting quark, 
yielding a quadratic sin((/:)) term in the angular depen- 



dence of not only the produced hadron but also the next 
remnant quark in the decay chain. Consequently, from 
this effect alone the distribution of each consecutive emit- 
ted hadron acquires another order in sin((/?). For a de- 
tailed investigation of this effect we refer to Matevosyan 
et al. J9], who employed a toy model to study the ef- 
fect of quark spin flip in the decay chain for unfavored 
fragmentation functions. 

The second source of the angular modulation arises 
from the spin flip probabilities themselves. This can be 
easily seen from the expressions for the spin flip proba- 
bility of Eq. ([3]), where = —p± and Px = p±cos{(p)^ 
Py = P± sm{ip). Thus both |aip and |a2p are modulated 
with a polynomial of the form a + 6sin^((^). Together 
these effects create an angular modulation of the hadrons 
produced at n'th order that is a polynomial in sin((/?) of 
order 1 + 3 x (n — 1). Then we can write the general form 
for the polarized fragmentation function for R hadron 
emissions as 

l+3(i?-l) 

D^/q,{z,Pl^)= Cn{z,Pl)sm^^. (4) 

Thus, in the case with only two hadron emissions, we 
need a fourth order polynomial to fully describe the an- 
gular modulation, as was shown in our fits in Fig. [i] 

Figure [5] shows the coefficients, c^, of the n'th power 
of sin((/:?), obtained from a fit to the fragmentation 
function for a u quark going to a 7r+ at z = 0.1 and to 
a 7r~ at 2: = 0.3, as a function of P^. In this case the 
Monte Carlo simulations involved six emitted hadrons 
(beyond which the results are stable). Here we do not 
show the constant term corresponding to the unpolar- 
ized fragmentation function, as it is much larger and has 
been extensively studied in our previous work [161 flTj . 
The linear term, ci, corresponds to the Collins function. 
The plots clearly show that the higher order terms fol- 
low the pattern set by the linear Collins term studied 
in Ref. [17 : they vanish for large values of P^ and os- 
cillate with an increasing magnitude with P^ ^0. In 
addition, the higher order terms are comparable in mag- 
nitude to the Collins terms in the small z region. It is 
worth noting, that the higher order terms for unfavored 
fragmentation remain significant compared to the regu- 
lar Collins term to substantially larger values of z than 
the favored fragmentations. 

We note that the hadrons produced further down the 
quark decay chain carry increasingly smaller momentum 
fractions and thus impact the distribution functions only 
in the low z region. This was shown explicitly for the in- 
tegrated, unpolarized fragmentation function in Ref. [15], 
where the plots in Fig. 5 demonstrate that the fragmen- 
tation functions hardly change after just three hadron 
emissions in the region z > 0.2. Thus the higher order 
modulations of the distributions appear mainly in the 
low z region, making it harder to isolate the higher order 
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FIG. 5. The fitted values for tfie polynomial coefficients in 
sm{(p), Cn, (a) for tt^ at z = 0.1 and (b) for tt" at z = 0.3, 
produced hy au quark for simulations with 6 emitted hadrons. 

polynomial terms in the MC data. Thus, a fourth or- 
der polynomial was sufficient to describe our MC results 
within the statistical errors with 6 hadron emissions. 

In this work we have described for the first time the 
higher order Collins modulation effects for the hadrons 
produced in the fragmentation of a transversely polar- 
ized quark within the quark-jet hadronization picture. 
We showed that these higher order polynomial terms in 
sin((/9) necessarily arise from two effects in each hadron 
emission step: the linear Collins modulation of the rem- 
nant quark number densities in the decay chain and the 
sm^{(f) modulation of the quark spin flip probabilities. In 
principle, this leads to an increase of the angular modula- 
tions of the resulting hadron distributions by three orders 
in sin((/9) with each hadron emission after the first one. 
We showed that with just a linear Collins modulation 
as an input in the elementary quark emission our Monte 
Carlo results require a fit with a polynomial of at least 
fourth order to achieve satisfactory values of Xdof • Earlier 
toy model tests in Ref. [19^ helped to establish this mech- 



anism of generating higher order modulations. Moreover, 
they showed that the experimental suggestion that the 
1/2 moments of the favored and unfavored Collins func- 
tions may have opposite signs, may be understood in the 
quark-jet picture if, on average, in the elementary hadron 
emission step the remnant quark's spin is antiparallel to 
that of the fragmenting quark. 

These new Collins modulation effects are only visible 
in the transverse momentum dependent hadron number 
densities and grow relative to the unpolarized and regu- 
lar Collins terms in the small z and region, where the 
multiple hadron emission effects are the largest. Never- 
theless, if observed experimentally, these effects will pro- 
vide a critical test of the quark hadronization process. 
Further, their inclusion in the analysis of high precision 
experimental data will allow us to improve the descrip- 
tion of the data and reduce the systematic errors in the 
extractions of the unpolarized and Collins fragmentation 
functions. 
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